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Abstract 

We discuss the stabilisation of the inverse cascade in the large scale instability of the 
Kolmogorov flow described by the complete Cahn-Hilliard equation with inclusion 
of (3 effect, large-scale friction and deformation radius. The friction and the fi values 
halting the inverse cascade at the various possible intermediate states are calculated 
by means of singular perturbation techniques and compared to the values resulting 
from numerical simulation of the complete Cahn-Hilliard equation. The excellent 
agreement validates the theory. Our main result is that the critical values of friction 
or (3 halting the inverse cascade scale exponentially as a function of the jet separation 
in the final flow, contrary to previous theories and phenomenological approach. 
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1 Introduction 



Inverse cascades are a common feature in the large-scale velocity and magnetic 
fields of geophysical, planetary and astrophysical two-dimensional flows. Their 
halting by spontaneous formation of zonal jets has been object of great interest 
and considerable work by many scientists, see among other Refs [1-5]. The 
phenomenon of an halted inverse cascade could play a role in the atmosphere 
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of Jupiter and other Jovian planets which exhibit jet streams of east-west and 
west-east circulation. 

Frisch et al. [6] showed that the inverse cascade of the large-scale nonlinear 
instability of the Kolmogorov flow described by the Cahn-Hilliard equation 
[7-9] may be stopped by the dispersive Rossby waves, i.e. by the so-called the 
/3-effect. We recall here that in the absence of any stabilizing effect the inverse 
cascade proceeds by visiting a family of metastable states with increasing 
scale until the final largest scale is reached [10]. In a later paper, Lcgras ct 
al. [11] have used singular perturbation techniques to calculate the range of 
the P values stopping the inverse cascade in one of the otherwise (if /3 — 0) 
metastable state. Their result was different from that obtained by the standard 
phenomenology based on dimensional arguments [1] which fails because it 
does not take into account the strong suppression of non linearities in the 
metastable states of the inverse cascade. 

In Refs. [6,11] the forcing maintaining the basic flow was chosen parallel to 
the planetary vorticity contours and there was no friction or advection. In a 
more realistic setup, Manfroi and Young [5] studied the stability of a forced 
meridian flow on a /5-plane pushing the fluid across the planetary vorticity 
contours, and including both friction and advection by a mean flow. They ob- 
tained a complete amplitude equation for the leading order perturbation, from 
which, with some formal modifications and with a slightly different interpre- 
tation of the parameters, the amplitude equation of Ref. [6] can be recovered 
(see Section. 2.1). Using this equation, but disregarding the dispersive contri- 
bution from /3-effect, they showed that random initial perturbations rapidly 
reorganize into a set of fast and narrow eastward jets separated by slower and 
broader westward jets, followed then by a much slower adjustment of the jets, 
involving gradual migration and merger. The stabilization discussed in Ref. [5] 
is only due to friction and not to dispersive effects. 

In this paper we present a systematic approach to the study of stabilization of 
the inverse cascade in the supercritical regime of the large-scale Kolmogorov 
fiow provided both by friction and P effect. The stabihzation by /3 effect al- 
ready discussed in [11] will be here presented in a greater detail and stabi- 
lization by friction will be discussed in the same mathematical approach as 
for the [3 case. The mathematical framework is based on the kink dynam- 
ics introduced by Kawasaki and Ohta [12] to describe the solutions to the 
Cahn-Hilliard equation; singular perturbation technique is used to study the 
stability of these solutions with respect to small perturbations due to fric- 
tion and /? effect. The perturbative calculations are performed analytically for 
large wavcnumbers and numerically for all cases. The results are compared 
together and with direct numerical stability and time-dependent solutions of 
the amplitude equation. 
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The paper is organized as follows: in Section 2 the Cahn-Hilliard equation 
in its complete form is presented; by complete in this context we mean that 
(3, friction, advection velocity and deformation terms have been added to the 
standard Cahn-Hilliard equation. The perturbation to the steady metastable 
solution of the Cahn-Hilliard equation by small (3 and a friction terms is pre- 
sented in Section 3. Section 4 discusses the stability of the steady metastable 
solutions and provides the main results of this work. Section 5 presents the 
techniques used to solve numerically the perturbation and the time-dependent 
problem, and compares the results of those calculations and the analytical re- 
sults. Section 6 offers a summary and conclusions. 

We do not give the detailed presentation of the mean advection effect, which 
introduces considerable additional algebraic complications, in order to focus 
here on the essential mechanisms that stabilizes the Cahn-Hilliard cascade. 
The main results in presence of mean advection are, however, given without 
demonstration in Appendix D. 



2 The complete Cahn Hilliard equation 

2.1 Basic equation 

Our starting point is the large-scale Kolmogorov flow in its slightly supercrit- 
ical regime, which is described by the Cahn-Hilliard equation [7-9] . We recall 
here that the basic Kolmogorov flow, u = (cosy,0), is maintained by a force 
f = i/(cos y, 0) against viscous dissipation. This flow exhibits a large-scale in- 
stability of the negative eddy viscosity type when the kinematic viscosity z/ is 
slightly below the critical value Uc = l/V^. 

Taking further into account /3 effect, friction r, external deformation radius 
l/S* Q and advection effect due to a non zero mean velocity 7 , we obtain 
by multi-scale techniques the following adimensional equation for the leading 
order large-scale perturbation : 



dt{l - S^d;^){v - 7) = \dlW\v) - A9> - pd-\v - 7) - rv, (1) 



where ^ denotes the integration in x defined for the family of functions with 
zero average over the interval [0,L]. The constants in (1) are 



The external deformation radius accounts for the inertial large-scale effect of a 
free surface in geophysical flows [13]. 
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(2) 
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and the potential W{v) is 



W{v) = - s'v' (3) 

Equation (1) was derived in [6] with S* = r = 7 = 0. In this derivation, 
the Kolmogorov basic flow is oriented in the zonal direction on the /3-plane 
and therefore the large amplitude flow develops in the meridional direction. 
Using a more realistic setting where the Kolmogorov basic flow is oriented in 
the meridional direction as an idealized baroclinic perturbation, and introduc- 
ing friction and mean advection velocity, Manfroi and Young [5] derived the 
following amphtude equation 

drA = -rA - (2 - 7*^) (9^2 A - 30^4 A 

+ 2Yd,{d,A^) + '^d,{d,A^) - Pd,-.A (4) 

where r and t] are the temporal and spatial variables. In their study, the 
last term on the right-hand side, which is the only dispersicve term in the 
equation, was set to zero. 0. By the change of variable dxA = w — 7 and the 
rescaling t = at, r] = bx, w = cv and 7* = 07 where a = 12s^A/(2 + 7*^)^, 
b = (6sV(2 + 7*2))i/2 and c = (3(2 + Y^)/{2T^)y/^, this equation is easily 
transformed into (1) up to the term in S which is only a trivial modification 
[13]. We see that the quadratic term in (4) disappears in this transformation. 
The 7 term represents the net advection velocity (see [5]) and is also the 
average value of v. 

In a recent work [14], it is claimed that the Kolmogorov instability exhibits 
a singular limit when /? <— O. This result is, however, established for a very 
special situation which does not arise here [15]. 

2.2 Kinks and antikinks 

The pure Cahn-Hilliard equation, is recovered from (1) by setting (3 = r = S = 
7 = 0. It admits a Lyapunov functional and is therefore integrable [8]. This 
property is preserved in the presence of friction but is lost in the presence 
of (3 [^. The solutions to the pure Cahn-Hilliard equation live essentially, 
albeit some initial transients, within a slow manifold of soliton-like solutions 

^ The coefficient /3 in the last term of the right hand side of (4) is the product of 
the planetary vorticity gradient per the small angle between the Kolmogorov flow 
direction and the planetary vorticity gradient. Therefore it can be either positive or 
negative unlike in the derivation of Ref. [6] where it is always positive 
^ A Lyapunov formulation for the approximated equation is recovered, however, in 
the limit of large /3 [6] . 
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Fig. 1. Kink-antikink annihilation in a numerical simulation of the pure 
Cahn-Hilliard equation with L=76.95 . 



with an alternation of plateaus v — ±r, separated by alternating positive 
and negative kinks, that we will call respectively kinks and antikinks [16] in 
the following. For large enough separation between adjacent kinks, the kink 
centered in a; = Xj is locally given by 

Mj{x) = ejM{x — Xj) = EjV tanh s{x — Xj) , (5) 

where ej = 1 for a kink and ej = —1 for an antikink [12]. This solution satisfies 
the equation 

-dlM, + W'{M,)^0. 



Within a x periodic domain of period L, the Cahn-Hilliard equation exhibits 
stationary metastablc solutions of period A = L/N with N pairs of alternating 
and equally spaced kinks and antikinks. These fixed points are unstable saddle 
points of the Lyapunov functional, except for = 1 which corresponds to an 
absolute stable minimum. The temporal evolution characterized by a growing 
total energy is a cascade of annihilations of kink-antikink pairs (see Fig. 1) 
leading eventually to the gravest mode iV = 1 [12]. It is shown in Appendix A 
that the local solution (5) is modified by terms of order exp(— sA) when the 
periodicity is taken into account. 
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3 Perturbation of stationary solution under the action of (3 and 
friction 

In order to study the modification of the stationary solutions to the pure 
Cahn-Hilliard equation under the effect of a small f3 or friction, we multiply 
both terms by a small control parameter, e, keeping (3 and r as 0(1). 

We put 7 = 0, leaving the case 7 7^ 0, which is technically much more involved, 
for Appendix D. We only notice here that 7 > breaks one symmetry and 
generates narrow intense westerlies and broad narrow easterlies [5] . 

It is convenient to integrate (1) twice in x to obtain 

_ 1^-2(1 - S%')dtv - e^d^h - e'-dfv = dlv - W\v) + h{t) , (6) 

where h{t) arises from integration in x. Then the perturbed solution is defined 
as U = v^'^'^ + ev^^^ + e^v^"^^ + 0(5^), where satisfies the stationary CH 
equation — (9^I;(°) + = 0. We will distinguish two cases for 5": (i) 5* = 

0, associated with synoptic and subsynoptic dynamics, and (ii) S = 0(1), 
associated with planetary motion (cf. [13]). When /3 7^ 0, we introduce the 
phase velocity c = eci + £^C2 + 0{e^) of the traveling framework in which v is 
stationary. 



3.1 Order 1 perturbation 



We first treat the case S — The first order perturbation v*^^) satisfies the 
linear equation 

t{v^'^) = g(°) (7) 

with 

Tig) = dig - <^ , 
ly; = W^"(^;(o)) , 

dxV^^^ belongs to the kernel of !F: this can be easily verified by multiplying (7) 
by dxv'^^^ and integrating within the domain. After integration over the spatial 
period, the solvability condition for (7) gives the first order contribution to 
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the phase velocity 



ci = ; ; — ■ (8) 



The symmetry of the stationary Cahn-Hilliard equation with respect to x 

reversal is inherited by v^'^\ The solution is antisymmetric with respect to 
kinks locations and symmetric with respect to the middles of the plateaus. 
The perturbation v^^^ is the sum of two parts 

where v^^^^ has the same symmetry as v^^^ and Vp^^^ has opposite symmetry. 

It is shown in Appendix A that, up to errors of 0{e~^^^'^) the basic solution 
can be approximated by a series of jumps locally described by (5). Over 
each interval [xj — A/4, Xj + A/4] centered on a kink in Xj (see (5)), we have 



from which the velocity ci is readily calculated using (8). After some algebra, 
we get 

where ( is the Riemann zcta function and A = ln^(l+e^^)(ia; = 0.150257 ■ • • . 
The phase velocity is directed to the left when f3 > and increases with the 
wavelength A. The fact that the error in (9) is exponentially small makes this 
expression very accurate even for not so large values of A as we shall check 
below. By straightforward algebra, one gets 



^ 2A ' ' I 3 4 ' ' 24 , 



-2aI"I 



(ixl -^a) +0(/3A2 + rA). (10) 



At distance from the kinks the ratio between the derivative and the potential 
term in is 0(1/A^). At leading order in 1/A: 

¥'^ = w^\4\^\-l^)(\^\-]A^ (11) 



24s2A' 'V ' 2 /V ' 4 
^^^^ = ^^\x\ (\x\ - 1-A) . (12) 



8s2A' 'V' ' 2 
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These expressions are valid at 0(A) distance from the kinks. It can be shown 
that at 0(1) distance from the kinks, v^^^^ is O(A^), while v^^^^ is 0(A). 

For S = 0(1), the leading contribution of friction is unchanged, but the effect 
of P is deeply affected. Q^^^ is modified as 

= i (c,d-'v^'^ - (/? + ci52)a-V°) - ra-%(°)) . 

Therefore, the phase speed is now 

where Ci is given by (9), that is 

<^^s--^,+j^,+0^j,j. (13) 

Unlike the infinite radius case, the phase speed varies weakly with A. Reporting 
(13) in Q^g^ and solving for v^''^\ we obtain 

at distance from the kinks. The correction to the stationary solution scales as 
A and is considerably reduced with respect to the case 5" = 0, for which it 
scales as A^. 



4 Stability 

4-1 Stability of the Cahn-Hilliard equation 

The stability of the solution = is easily obtained by linearizing (1) in the 
Fourier domain. When r = 0, the solution is unstable to all Fourier modes 
with wavcnumbcr Q < k < km = V^s = ^2/3. This result docs not depend on 
the values of f3 and S. When r > the modes near k = and km are stabilized 
and the instability band in /c-space shrinks as r grows. The solution v = is 
stable for r > tq = Xs'^ = (3\/2)~^. The vicinity of this value has been studied 
in [5]. We are here interested in the limit of small r. 

The stability of the non-zero stationary solutions to the Cahn-Hilliard equa- 
tion can be studied using the equation for kink motion derived in Appendix B. 
For an arbitrary perturbation, fast transients dissipate rapidly, leaving only 
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after a short time the part of the perturbation that projects onto kink dis- 
placement. The jth kink being displaced by dxj, the perturbation to v^'^^ is 

2N-1 

Sv^- J2 d.MeSxe. (15) 



Then using (B.14), the equation for the displacements is 

A-p2 2N-1 / 2 1 \ 

= 64:S^r^e-'^[2Sxj - Sxj+i - Sxj-i] - 2re//i . (16) 



It is convenient to use Fourier components defined as 

2N-1 

SXj = Yl V'me^^^ , 

m=0 

with m e [0, 2A'" — 1] and 4'2N-m — V'm- The equation for is obtained 
by multiplying (16) by and summing over the j from to 2N — 1. 

Taking into account the regular alternation of kinks and antikinks separated by 
intervals of length L/2N in the basic solution and using the Fourier transform 
of the Green function given by (C.2), one obtains after some algebra: 



■(—^^—- --)lP^= 1285^6-^^(1 - 
.\1+C0s9m Sj 



N 

- 4:—Sh 5{N - m) , (17) 



with 9m = nm/N. The leading order form of (17) was given by Kawasaki and 
Ohta [12] with a factor 2 error (see Appendix B). 

The first term in the right hand side of (17) is destabilizing the stationary 
solution with the eigenvalue 

1285^Ae--^ . , ( 2(l + cosgJ V^ 
(To = -j^ sm ^^11 -j^ 1 . (18) 

This instability is responsible of the inverse cascade in the CH equation. 
Each value of m is associated with a real eigenvalue of T and a dimension 
2 eigenspace. It turns out that an appropriate basis of this eigenspace is pro- 
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(a) (b) 




Fig. 2. Numerical calculation of Va and Vi, by discretization of the eigenvalue problem 
for (of Section 5.1) with 128 points for L = 76.953 and m = A. (a) thin solid: 
v^^'^; solid: dash: vi,; (b) thin solid: dxv'^'^'^; solid: Va/dxV^^\ dash: Vb/dxV^^^. The 
scale is arbitrary for Va and Vb- Even if the plateaus in v^^^ are very short for the 
chosen values of L and m, the ratios in (b) show the staircase structure of Va and 
Vb over the kinks as given in (19, 20) except where dxV^^^ vanishes in the middle of 
the plateaus. 

vided by the couple of orthogonal vectors 

2N-1 

M^) = Yl i-^y cos jemdxM{x - Xj) , (19) 

j=0 
2N-1 

Vb(x) = (-1)' smjOmOxMix - Xj) , (20) 

j=0 

which are here given up to an error 0(e~*^/^). These expressions agree very 
well with the numerical solution shown in figure 2. 

The 6h contribution vanishes but on the mode m = N. The solution v{x,t) 
must average to zero within the periodic interval [0, L] . In terms of kink mo- 
tion, this imposes the constraint 

and thus, ipN — 0. The presence of 5h{t) in (17) is required to impose this 
condition. 



4-2 Stability of (3-CH equation with infinite radius of deformation 
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4-2.1 Formulation of the problem 

The equation governing the perturbation dv to v^^^ can be written conveniently 
for (fi — d~^5v 

dtif^C^, (21) 

with 

C = -XdM - W"{v))d, - cd, - el3d-^ - er . (22) 

The reason of using ip instead of 5v is that L is auto-adjoint while the corre- 
sponding operator for 5v is not. 

Unlike the case /5 = 0, the slow component perturbation to the stationary 
solution of the /3-CH equation does not reduce to the simple motion of kinks. 
One has also to take into account the dispersive effect of the term modifying 
the shape of the slow modes and contributing to the stability. Therefore, we 
expand C as 

£ = £0 + + + O (£') , (23) 

with 

C^^-\dM-Ul^)d., (24) 
>Ci = A9,«%(^)9,) + cA - Pd-' - r , (25) 

£2 = Aa,(iy^V2) + + C2d, . (26) 

The eigenvalues of (21) are perturbations of the eigenvalues of (16). For a 
given m ^ N, we obtain 

a ^ ao + sai + ie/ii + ea2 + ieii^ + Oit") . (27) 

The functions = dx^^a and 995 = d^^vi, are orthogonal eigenmodcs of Co 
and, it turns out, an appropriate Jordan basis for the perturbation problem. 
They are respectively modified as ipa + eipai + £'^(Pa2 + 0{e'^) and ipb + ecpbi + 
s'^<Pb2 + 0{e^). The hierarchy of hnear problems is 

Co'fa = cro(pa , (28) 

CoiPal + Ciipa = (JQiPal + aiipa-f^lfb , (29) 
jO-0<Pa2 + Ciifial + jO,2(Pa = CF0^a2 + CFi^Pal + a2^a - l^l<Pbl - l^2^b ■ (30) 

and similar equations for ipi^- 
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4-2.2 Stabilization by friction at first order 

The first order corrections of tfie eigenvalue a are obtained as solvability con- 
ditions of (29) by 

< (Pa,Ci(pa> ^(Ti <iPa,Va> , (31) 

< ipb, CliPa > = -A*i < Vb, > , (32) 

where < f,g >= Jq f{x)g{x)dx. In < (pa,Ci(pa >, the contributions from P 
immediately vanish by integration. We have, after integration by part 

< ^a, J^l^a >- -A < Va^ > +r < Va, dfva > ■ (33) 

The first contribution to the right hand side of (33) can be reduced to an 
integral over a single kink by summing the trigonometric factors arising from 



< vj, >=jj < {d^MjYdx . (34) 

Here j labels an arbitrary kink and the integral bounds do not need to be 
specified since d^Mj decays exponentially on both side. This contribution is 
further transformed using 

j W^'v;^^\d^Mjfdx ^jj d-'v^'^^d^Mjdx, 

which is valid up to exponentially small errors. Finally the contribution is 
reduced to a non local integral by part, 

1 /• 1 .n^„ - , , 1 ^ .n^„ , /n^ , l„o 2r2 



As 



Combining this with (C.8) and (C.2), we obtain 



(ji = -rsin^ ]^9m 



+ £ (l + cos^ -Om) (l - ^ cos^ ^Orn) + 0{e''^/^) . (35) 



In a similar way, we have 

+ Ci < Va, d~^Vb > -P < Va, d~^Vb > ■ (36) 

The first term in the right hand side of (36) vanishes after the trigonometric 
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summation. Using (C.8-C.9) and (C.1-C.3), we obtain 

x(y(l + ^^)-^)" + 0(e-^/^), (37) 
or, retaining only the first three orders of the expansion, 

'"'-^y 12(1 + 3s(l + ^2)2 + 3As2 (1 + ^^)3 j^UV' 

with t = tan7rm/2A''. 

It is interesting to notice that the nonhnear contribution < Vo^jUq'v^^^ > is 
destabihzing the stationary solution. However, the direct hnear damping by 
friction remains larger and the total effect of friction is always stabilizing. 
Therefore, the m-mode perturbation to the stationary solution is stabilized 
by friction for 

^.^e~^^ o, , irm , 

r > Tc = 512^— s^A cos^ — , (39) 
A IN ^ ' 

at leading order. 

The corresponding result when 7 7^ is given in (D.6). The mean advection 
decreases the value of friction necessary to stabilize a given wavenumber. 



J!f..2.3 Stabilization by (3-effect at second order 

At second order in (72 is solution of the solvability condition. Here wc set 
r to zero for simplification as stabilization by r is already obtained at first 
order. 

< Va, ^Wal > + <Va, ^2<Pa > C^a < (fa, ifa > -j^l < Va, Vbl > ■ (40) 

The second term on the left hand side of (40) expands as 

< Va, C2Va > = -A < va', ^V^^ > -^A < <^(IJ(^))^ > . (41) 

The second term is 0(/3^A"^) while the first term is 0(/9^A"^) and dominates 
at leading order. This term can be transformed in the same way as above for 
< (Pa,J^i(Pa >, leading to 
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where Q^^^ is the right hand side for the second order version of (7). After 
dropping out all contributions that vanish owing to the symmetries, we obtain 

+ 2X / <"i'^^'^)"9lM,dx + 0{f3'A') . (42) 

The second term on the right hand side of (42) is negligible relatively to the 
first. This latter can be calculated using 

AJ ^ 2AA,-iA ^ ^ 2AA,-iA 

and (11). We obtain 



The other contributions in (40) involve ipai and ipbi. These quantities can be 
approximated in the same way as v^^^ at distance from the kinks. We have 

is'Xcpai = -A9-^«'tj(i)^;„) - ci9jVa + - f^lO^'^b , (44) 

and a similar expression for (pbi. Then, after some algebra and dropping the 
terms which do not contribute to the leading order, we have 

- 2/?Ci < Va, d-\a > -flip < Va, d'^Vb > +filCi < Va, d'^Vb > +0(/3^A^) . 



Using the results of Appendix C, we obtain 

92,160 s^X 



and 



A3 



Finally 

a - ^!^^!^i±^ + Or/5^A3l (45) 
- "69,120 .^A {l + t-y + 0(/3 A ) . (45) 

The contribution from is a correction to the first order stabilization obtained 
in (7i. The (5 term does not appear at first order and is stabilizing in (45). 
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Though the effect is small, it increases algebraically with A while the non 
linear coupling of kinks decreases exponentially in (18). Therefore, if r = 0, 
stabilization of the m-mode perturbation to the stationary solution is obtained 
at leading order for 

/5 > /5e = (35,389,440 -^s^A^^— ^ j . (46) 
The condition is the most restrictive for m = 1, that is t — tauTr/ (2iV). 



4-3 Stability of (3-CH equation with finite radius of deformation 



In this case Ci and C2 are modified as 

Ci = Aa,<V^)9,. + + c,S^)d-^ - r , (47) 

A = A9,«V2) + ^[/o'V^)>, + C2(9, - S^d-') . (48) 



At first order in e, the imaginary part of the eigenvalue is 



4(3t{2 + 3t') ^ J P 



(49) 



At second order in £, (72 is still given by (40). The contribution < (pa-ii^2^a > 
is now 

After integration by part and using (14), we obtain 

Similarly, we have 

< > = -^^^ (1 + + 15^^) + O (xj ' 
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Finally, we obtain a-i as 



The stability crossover for /3 is now 

/Se = (l5,360 s^5^A2^i-^^±^^) . (51) 

Therefore, the stabilizing effect of /3 is much reduced compared to that with 
infinite radius of deformation. 

All the perturbative calculations of Sections 3 and 4 have been checked with 
Mathematica. 



5 Numerical approach 



The analytic results established in Section 4 are valid in the double limit of 
small e and large A. These asymptotic results are complemented and com- 
pared with three types of numerical calculations: (i) numerical solution of the 
perturbative problem for several values of A, (ii) direct numerical simulation 
of the Cahn-Hilliard equation in the Fourier space and (iii) direct stability 
calculation. 



5.1 Numerical solution of the perturbative problem 



We relax here the hypothesis on the large value of A by solving numerically the 
perturbative equations for v^^\ v^'^\ cpai and cpbi, and numerically evaluating 
the solvability conditions. 

The calculation is performed according to the following algorithm 

(1) The basic solution v^^^ is defined by the approximate form given in ap- 
pendix A. 

(2) The inverse derivatives d~^v^^\ and d~^v^^^ are calculated by 
Fourier transform and tabulated to obtain Q^^\ 

(3) Ci is calculated by discrete evaluation of (8). 

(4) (7) is discretized as a tridiagonal problem and solved for vp^^^ and v^^^\ 

(5) Q(^) is built in the same way as Q*-"-* from the inverse derivatives of vp^^^ 
and TV . 



16 



(6) v^"^^ is calculated by inverting a tridiagonal discretized problem. 

(7) The eigenvectors Va and are defined using (19) and (20). 

(8) ipa and ipb are calculated by Fourier transform and tabulated. 

(9) (Ti and /^i are calculated according to (31), (32) and (25). 

(10) (pai and (pbi are obtained using (29) and solving a tridiagonal discretized 
problem. 

(11) (72 is obtained from (40), (25) and ((26). 

This algorithm admits 0{e~^^^^) errors, but all the steps generating errors of 
algebraic order in the asymptotic expansion are here solved numerically. The 
implementation has been done as a Mathematica notebook available from the 
authors. The number of grid points and Fourier modes has been adjusted as 
a function of L and N in order to provide at least 3 digits of accuracy in the 
results. The symmetries have been exploited to distinguish the solutions and 
reduce the number of points. 



5.2 Numerical solution of the complete Cahn-Hilliard equation 



We have done numerical simulations of the complete Cahn-Hilliard equation 
(1) for the case 7 = S* = 0. For practical convenience, the spatial period has 
been kept fixed to 27r by rescaling x as x ^ px. The complete Cahn-Hilliard 
equation then reads 

drv = ^9lv^ - ^9lv - ^d^v - pPd'^ - rv . (52) 

Since Fourier modes are discretized by the periodicity condition, the number 
n of unstable modes ior v = is the integer part of p(2/3)^/^. 



5.2.1 Time- dependent simulations 

The simulations are performed using a standard semi-spectral method where 
the number of retained real Fourier modes is 256. The collocation grid in the 
spatial domain has 512 points in order to fully remove the aliasing due to the 
cubic term in (52). We have checked that using higher resolution does not 
modify the results within the explored parameter range. The temporal inte- 
gration is performed with an Adams-Bashforth second-order scheme. Initial 
conditions are a random white noise in the spatial domain. We have made 
a large number of runs by varying the initial conditions (changing the seed 
of the pseudo-random number generator and the amplitude of the noise) , the 
values of r and n. 

Section 4 shows the existence of multiple stable solutions induced by P and 
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friction but does not provide indications about the attraction basins of these 
solutions. In the inverse cascade of the standard Cahn-Hilliard equation, the 
interaction of a pair of neighbor kink and antikink scales as exp(— sAo;), where 
Ax is the distance between the kink and the antikink (see Appendix B). 
When friction r is just above the critical value rc{N) stabihzing the solution 
with N pairs, wc may conjecture that the stabilizing effect is of the order 
0((r — rc{N))6x), where Sx is the departure of a kink from its equilibrium 
position (this is clear from the shape of the eigenmodes Va and (cf. ection 
4.1)). This effect, however, cannot extend very far in Sx as the attraction to the 
neighbor antikink grows as exp{s5x). Moreover, the time-dependent solutions 
do not need to pass in the vicinity of the A'"-pair fixed point during the cascade 
of kink-antikink annihilations. Therefore, we expect that the fraction of the 
solutions halting on the A^-pair stable state will be small in the vicinity of the 
critical Tc and will be significant only when the stabilizing effect is felt over a 
distance of order A. Once this is obtained, very few solutions should jump to 
lower N states. 

Figure 3 shows how the inverse cascade evolves as a function of r for the same 
initial conditions, with the corresponding steady states shown in figure 4. The 
final state wavenumber increases with r and stays bounded by the value of 
the most unstable eigenmode of the Kolmogorov flow k = (2/3y^^n which is 
the unique mode excited when r approaches Tq. The total energy calculated 
as the sum of the various energies E(k) of the single modes k, decreases as the 
friction increases. 

Figure 5 compares the halting of the inverse cascade by friction and by /3 effect 
for the same initial conditions. It is apparent from flgure 3 that friction simply 
halts the cascade by stabilizing one of the intermediate steps. The effect of the 
/3-term is more complex: oscillatory transients are excited and, paradoxically, 
the cascade is accelerated. In Fig. 5 the transitions to A" = 3 and to A" = 2 
occur much earlier than in the absence of beta. Other examples can be found 
in [17,11]. In the final state, the f3 effect breaks the symmetry between the 
kinks and the antikinks which is preserved by friction. 



5.2.2 Numerical stability calculations 

In order to compare the results from the analytic and numerical perturbation 
approaches to the direct simulation of the complete Cahn-Hilliard equation, 
we have developed a direct analysis of stability by the same technique as for 
the time-dependent simulations. The core of this analysis is to calculate the 
Jacobian matrix of the right-hand side of (52) linearized around a given state v, 
with respect to each of the Fourier component. This is done by differentiating 
(52) and calculating the columns of the Jacobian matrix by the semi- spectral 
method applied to the differentiated equations for each Fourier component. 
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Fig. 3. Temporal evolution of the energies E{k) of the Fourier modes and of their 
sum Etot- The three shown cases, are for r = 0, r = 310^ and r = 10^, with 
/? = 0, n = 10 and the same realization of white noise as initial condition. For 
r = the inverse cascade is complete to iV = 1; for the other values the cascade 
stops respectively on A'^ = 2 and on A'^ = 3. Note that Efot is constant between two 
annihilation events. Increasing further r enables to stop the cascade on larger N 
configurations (not shown). 
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Fig. 4. The corresponding asymptotic velocity profile for the three cases presented 
in Fig. 3. Note that the kinks-antinkinks pairs arc equidistant, unlike the iV = 2 or 
N = 3 configurations of Fig. 1 which are only metastable. 

The stability calculation is performed as follows. First an estimate of the 

stationary solution based on Appendix A is refined by a Ncwton-Raphson 
algorithm which usually converges within a few steps. The degeneracy due to 
the x-invariance of the complete Cahn-Hilliard equation is removed by setting 
to zero the imaginary part of the dominating Fourier mode and removing the 



19 




Fig. 5. Comparison between two temporal evolution of £?tot and E(k), for the same 
initial conditions and for very similar final solutions with same N and final energy 
Etot ~ 0.5535. Top: r = 10"^ and /3 = 0. Bottom: r = and /? = 9.7 10"^ 

corresponding row and line from the Jacobian matrix. In the case /? 7^ 0, 
the solution is stationary in a frame traveling at a velocity c. This phase 
velocity is treated as an additional unknown increasing the dimension of the 
Jacobian matrix by one. The stability of this numerical solution is then found 
by finding the eigenvalues and the eigenmodes of the Jacobian matrix using a 
standard QR algorithm from LINPACK. In this procedure, we find both the 
slow components associated to kink dynamics and the highly damped modes 
associated with fast relaxing transients. The conditioning of the eigenvalue 
problem gets very bad as A increases as a result of the large separation between 
slow and fast eigenvalues, thus limiting the parameter range for the numerical 
calculation of stabihty. There is enough overlap, though, with the vahdity 
domain of the perturbative theory to provide detailed comparison. The number 
of Fourier modes used in this analysis has been 128, 256 or 384, depending on 
the values of L and N. 



5.3 Comparison of stability results 



Table 1 compares the critical values of friction and (3 estimated from the 
analytic perturbative expansion, the numerical solution to the perturbation 
problem and the numerical stability analysis. In the numerical stability anal- 
ysis the value of the parameter is adjusted by dichotomy from two values 
bracketing the transition. There is an excellent agreement between the three 
values of critical friction when A is large. At the largest values, however, the 
QR algorithm fails to converge and no results are obtained for the numerical 
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stability. When A is not large, that is when the kinks are not distant enough 
to neglect the contribution of (9(e^*^/^), significant discrepancies occur be- 
tween the different estimates. We see from the table that this occurs when 
> 3 10-2. 

There is also an excellent agreement between (3^*^^^ and for the same 

range of A values as for Vc- The analytical prediction (46), however, provides 
only an order of magnitude and is wrong by at least a factor two when the 
two other quantities agree by four digits. The reason is that the error in (46) 
depends algebraically on A unlike the error in (39) where the error exhibits 
an exponential dependence. Very large values of A are required to make (46): 
for n = 200 and = 5, we obtain /^p^'* = 2.24 lO^^^ ^nd = 2.35 lO^^^^ 
while for A^ = 2 we obtain /^p^""* = 1.842 IQ-^^^ and = 1.808 lO^^^^ This 
difficulty with P is entirely due to the need to solve completely the perturbation 
at order 1. The phase speed Ci and the frequency /ii, which are obtained as 
solvability conditions at order one, are known with the same accuracy as Tc, 
as can be checked in Table 2. 



5.4 Comparison of stability properties and time- dependent solutions 

In order to assess the distribution of final states among the multiple stable 
steady states we have performed ensemble simulations for a number of values 

of r and j3. Each numerical integration of the time-dependent numerical model 
described in Section 5.2.1 is characterized by n, r, j3 and the initial condition. 
We choose for this latter a white noise with amplitude A in the spatial domain. 
For each value of the parameters and for two values of the amplitude, A — 
0.1 and A — 1., we performed an ensemble of 100 independent simulations 
by varying the seed of the random number generator. After some time, all 
simulations converge to a final stationary or uniformly traveling state (if /9 = 
or /5 7^ 0, respectively). The non convergent cases are due to fast transients 
leading to nonlinear numerical instabilities. 

Table 3 shows the dependence on r when /9 = and n — 20. The distribution 
of final states agree qualitatively with the analysis presented in Section 5.2.1. 
Stationary states associated to given value of A^ are only reached for r larger 
than the critical value rc{N). The proportion of solutions reaching these sta- 
tionary states is about 15% for A — 0.1 and 10% for A — 1. when r/rc{N) fa 3. 
This proportion grows rapidly as r increases further while the number of states 
reaching solutions with smaller A^ falls dramatically. In practice, it is difficult 
to find values of r where more than 3 steady states are obtained. It is also 
visible that larger amplitude of the initial conditions favor smaller final A^ and 
more dispersion of the final states. 
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Table 4 shows the corresponding results as a function of j3 when r = 0. They 
are qualitatively similar to the precedings although up to 7 different final 
states are now observed for /5 = 1. 



6 Summary and conclusion 

Wc have investigated the stabilization induced by friction and j3 effect in the 
inverse cascade of the large-scale instability of the Kolmogorov flow. This prob- 
lem has been treated by solving the unidimensional complete Cahn-Hilliard 
equation using both numerical simulations and perturbation techniques 

In the standard Cahn-Hilliard equation, the kinks and antikinks are coupled 
by interactions decreasing exponentially as function of their separation. The 
number of kinks decreases with time and their average separation increases 
as a result of the inverse cascade. Therefore the Cahn-HiUiard coupling also 
decreases with time. 

The basic effect of friction is to damp uniformly the motion of the kinks and 
antikinks. Nonlinear effects due to the deformation of the kinks tend to re- 
duce this damping but cannot invert it. When the separation is large enough, 
the damping overcomes the destabilizing Cahn-Hilliard coupling halting the 
inverse cascade before it reaches the gravest mode. The dependence of the 
critical friction Vc upon the kink separation A/2 scales as Tc ~ (e~*'^/A). The 
perturbative approach yields a very accurate analytical expression of Tc be- 
cause the leading contribution to TcB*^, which is algebraic in A can be obtained 
exactly leaving out only terms which are exponentially small in A. 

Stabilization by /3-effect is more complex. It does not contribute to the lin- 
earization of Cahn-Hilliard equation around a steady state and appears only 
at the second order of the perturbative expansion in (3. As the first-order step 
of the perturbation expansion is only solved for the leading contribution in a 
1/A expansion, the analytic expression of the critical j3c is fairly inaccurate 
when compared to numerical solution of the perturbative problem or to di- 
rect stabihty calculations. Its scaling 13c ~ (e~*^/A^)^/^ is, however, correctly 
predicted. 

The presence of a finite radius of deformation 1/S , which slows down the 
fastest Rossby waves, is to provide less efficient stabilization by the (3 effect. 
The critical (3^ then scales as (3^ ~ (e'^*~*^/A)^/^. 

In the presence of mean advection, as in Ref. [5], the eastward jets are narrow 
and strong while the westward jets are broad and slow. The critical value of 
friction is reduced (see (D.6)) and scales as Tc ~ {e'^^^^~^/^'' / h) 
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Stabilization is demonstrated near the modified steady states of the Cahn- 
Hilhard equation. The attraction basin of these steady states depends on the 
stabihty of time-dependent solutions and has been investigated numerically. 
The results suggest a fairly simple pattern where, for most values of the pa- 
rameters, the phase-space is filled by the attraction basins of only 2 or 3 sta- 
tionary solutions. The boundary of these basins might be very complicated, 
even fractal. 

We have observed in the numerical simulations of the pure Cahn-Hilliard equa- 
tion that the inverse cascade does not always begin by the most unstable state 

= km] however it is in principle always possible to stop the cascade at such 
a state by enhancing the friction. The same is not true for j3, which is not 
always able to stop the cascade at the scale corresponding to the most unsta- 
ble state. As already noticed, P dispersive effect is more complex than friction 
effect. While the inverse cascade prior to the halting by friction does not differ 
from the pure Cahn-Hilliard case, the f3 effect is paradoxically accelerating 
the cascade before halting. We speculate that this is due to the propagation 
of fast Rossby waves superimposed to the kinks and increasing their coupling. 

Our result for the critical friction Tc differs from that given in Ref. [5] where the 
authors found Vc ~ A~^. Their reasoning was based on varying and minimizing 
the Lyapunov functional with respect to A. Our interpretation is that this is 
questionablesince the Lyapunov functional can only be minimized with respect 
to the solution, not with respect to the parameters. 

The scaling of the critical f3c also differs from the scaling (3^ ~ A~^ which 
would arise from standard phenomenology [1] by balancing the nonlinear and 
the dispersive term in (1). The reason lies in the suppression of nonlinearities 
in the slow manifold for solutions of the complete Cahn-Hilliard equation once 
the initial transients have been dissipated. In more realistic two-dimensional or 
quasigeostrophic flows, the presence of strong dominating jets and/or coherent 
eddies is similarly inducing a reduction of nonlinearities with respect to a plain 
dimensional estimate. This is why the prediction of a k"^ energy spectra [18] 
is usually not observed in forced two-dimensional flows [19] with the possible 
exceptions of the smallest quasi-passive scales of the motion. 
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A Approximate solution of the Cahn-Hilliard equation 



The stationary form of the Cahn-HiUiard equation d^v — U'{v) — can be 
recasted as 

d.U=^d,{d,vf, (A.l) 

from which the solution is obtained by quadrature. By integrating (A.l) once, 
we obtain 

U^^{d,vf-C, (A.2) 



where C is a constant which determines the value of dxV on the kinks and the 
periodicity of the solution. For the single-kink solution (5), we have C = |s^r^ 
and the asymptotic value of U is — C. Let us now assume 

where /i is assumed a small perturbation and try a solution under the form 

V — rtanhsx + /iv . (A. 3) 

By replacing in (A.2) and expanding, we obtain, at first order in /i, 

2dxV — —Asv tanh sx — Ts cosh^ sx . (A.4) 
This equation can be solved as 

V — —-^ f (3 + cosh^ sx) tanh sx H ^^^^ — ^ , (A. 5) 

16 \ cosh sx J 

where the condition v{0) = has been used. We can relate /j, to the period of 
the solution by using dxV = for a; = A/4, leading to 



/X I -sF ^tanh^ sa;(3 + 2 cosh^ sx) + 



3s tanh sx 
cosh^ sx ^ 

1„ sT 



1 1 si 
--Fscosh^sx ^ + s 2 — = 0- (^-6) 

2 J cosh sx 

When A is large the main contribution is 

fx = 64e-^^ . (A.7) 



Comparison with numerical solutions of (1) shows that (A. 3) with (A. 5) 
and (A.7) approximates periodic stationary stationary to less than 0.2% for 
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A > 10. One can build a solution over the whole domain by using (A. 3) over 
contiguous intervals containing a kink matched at mid-distance between ad- 
jacent kinks. The approximate solution is continuous and the discontinuity of 
its derivative at matching points is 0{exp{—s\xAK — ^kD) where \xak — xk\ 
is the distance between adjacent kinks. 



B Kink motion in the Cahn-Hilliard equation 

This Section is adapted from Kawasaki & Ohta [12] and corrects one error 
found in this paper. 

We assume that the solution is a combination of kinks which are individually 
described by (5). As d^Mi^x) and W'{Mi) decay rapidly away from x = Xi, 
the solution in the vicinity of the j-th kink is the sum of Mj{x) and of small 
contributions from adjacent kinks which are the small deviations from their 
asymptotic values. We write 

v{x,t) = Mj{x) +Vj{x,t) , (B.l) 

where Vj is small in the vicinity of the j-th kink, but takes finite value at 
distance. A valid expression in the vicinity of neighbor kinks is 

Vji^, t) = T.iMe{x) - M,(+cx))) + J2iMi{x) - M,{-oo)) , (B.2) 

where Mi{+oo) and Mi^—oo) are the asymptotic values at infinity for the basic 
kink profile. Here we assume that kinks and antikinks alternate(i.e. ejej+i — 
— 1), and that they are numbered from to 2N — 1 within the periodic interval 
[0,L]. The time dependence is entirely contained within the positions of the 
kinks {xj{t)}. 

The temporal evolution of the solution is then given by: 

2N-1 

dtv{x, = - E Xi{t)d,Mi{x) , (B.3) 
e=o 

where v{x, t) is governed by 

-^d;^dtV^dlv-W'{v) + hit). (B.4) 
A 

The function h{t) arises from the integration in x. The other terms arising 
from the integration vanish owing to the periodicity in x. 
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In order to estimate the motion of the j-th kink, we use Mj as a test function 
by multiplying (B.4) and integrating over the domain. Then we obtain: 



1 2N-1 L fL 

/ T 'E xAMjd-^Medx = / ((^.v - W'{v))d,Mjdx + / hdMdx . 
A B C 

(B.5) 



2N-1 L r-L 

A=--yxi / d^Mj(x)g2(x-x')d^Me(x')dxdx' , (B.6) 
X Jo Jo 



Contribution A can be written as 

e=o 

where Q2 is the Green function solution of 

-dlg^ix) = 5{x) . 

dxMj and dxMg are two well separated functions which contribute to the 
integral in (B.6) respectively in the close vicinity of xj and x^. By expanding 

^2('^ — x') near Q2{xj — Xf) and summing local contributions using J{x — 
XjfdxM^dx = (-l)%2r/6s2 and J J \x - x'ld^MjdxMedxdx' = AT^s, we 
obtain 

4r2 2N-1 / 2 -1 \ 

>1 = -— E - xe) + (-I)^-'y^ - . (B.7) 

Using the expression for Q2 within the interval [0, L], one obtains 
, 2Lr''^\ f{{x,-x,)[L]f {x,-xe)[L] 



'Y.\-iy-% _ {x^ 

e=o \ 



-iy-%-e] . (B. 



1 1 



6 24L2s2 4Ls 

Contribution B is expanded using 

W'{v) = W'{Mj) + W"{M^)v^ + W'j^^{M^,Vj) . (B.9) 

Like Vj, W'^]^ is small in the vicinity of the j-th kink but finite at distance. 
We have 



B^- W'^LdxMjdx + {dlMj - W'{Mj))dxMjdx 

+ [\dlvj-W"{M,)vj)dxMjdx. (B.IO) 
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The second integral in (B.IO) vanishes and the third one vanishes also after 
integration of its first term by part. Therefore, we are left with 



B = - j^W'^L^xM^dx. (B.ll) 
Using (3) we find 

In the vicinity of of the order of the tails of Mj+i{xj) — Mj+i{—(yD) and 

Mj_i{xj) — Mj_i{oo), that is 0(max(e~^*(^^+i'~^^\ Q-2s{xj-xj-i)-^_ ^^le vicinity 
oixj+i and vj is 0(1) while d^Mj is o(e-2<^.+i-^.)) and o(e-2<^-^i-i)). 
Therefore, the two main contributions to B in (B.ll) arise from the vicinities 
of Xj+i and Xj-i. In the vicinity of xj+i we use 

Vj = Mj+i + Sj+iF 

so that 

{;|(3Mj + Vj) = ejr^(2 - tanhs(a; - Xj+i))(l + tanhs(x - Xj+i))^ . 
Using also 

and replacing in (B.ll) with similar contributions from the vicinity of Xj^i, 
we obtain 

B = -32s2r2(e-2<^^+i-^^) - e-2^(^^-^^-i)) . (B.12) 
Finally, contribution C gives 

C = 2ejrh{t) . (B.13) 



Summarizing the results, one gets 

A \ L 

= 32s^r2(e-2^(^^+i-^^) - e-2^(^^-^^-i)) - 2ejVh{t) . (B.14) 

Eq. (B.14) shows that two neighbor kink and antikink attracts themselves 
A stationary solution is obtained when B vanishes for all values of j. This 
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condition is satisfied if the kinks and antikinks are equispaced over the interval 
[0,L]. Then, h{t) = 0. 



Eq. (B.14) can be used to determine the motion of the kinks far from the equi- 
hbrium, under the condition that the kinks and antikinks remain far enough 
to satisfy the approximations of the above calculation. 

The calculation of Kawasaki & Ohta [12] slightly differs from our own and 
is hmited to the leading order. They fail to take into account the exponential 
variation of vj near xj+i and xj^i. Therefore, their result for the leading order 
of B contains an error, being too small by a factor 2. 



C Green function in the periodic domain and calculations of cou- 
pling coefflcient 

Within the periodic domain [0,L], the 5 function is made periodic by adding 
a constant value —1/L everywhere but at the origin. One can also use 




The solution to 



92' 



IS 




where 



41] -I 

i((,x[l]r-a:[l] + i), 

-i-(6(x[l]f-15(x[l])^ + 10(x[l]f-x[l]). 



96{x) 



where x[l] means x modulo 1. 
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The calculation of the perturbed motion requires to calculate the Fourier trans- 
form 



27V-1 . - A 



which can be written as a function of 



2N-\ / n \V 



j=0 



For N, Sm{p) can be calculated using the following relations 



Wc get 

Sm(0) = 0, 
Sm(l) = 



1 + a ' 
1 

s^(2) = -,r— + 



l + a A^(l + a)2' 

1 3a 3a(l — a) 

gm(3) = -^r—+ .^rn , + 



l + a 2Ar(l + a)2 4A^2(l + a; 



3 



1 2a 3a(l — a) a(l— 4a + a2) 
s (4) — 1 \ 1^ 1 I 1^ i 

""^ ' l + a N{l + ay 2m{l + af 2^(1 + a)^ ' 
5a(l - lla + lla^ - a^) 



Sm(5) = h 

Sm(6) = h 



3a(l - 26a + 66a^ - 26a^ + a^) 
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where a — exp{—i9jn) with 9^ — nm/N . We also have 

sjv(O) = 2iV, 

sn{2)^^{2N-1){AN-1), 
«iv(3) = ^(27V-ir, 

'^^^^ ^ 24^^^^ ~ ^^^^^ ~ - 67V - 1) , 

^^(5) = ^(27V - 1)2(87V2 - 4iV - 1) , 

^^(6) = l34iiV^^2Ar - 1)(47V - 1)(487V^ - 487V3 + 67V + 1) . 

Using these relations, the Fourier transforms are readily calculated. For m ^ 
N, we have 

Gi{m) = s^{l) - ^Sm{0) + ^ = —it , (C.l) 
^2(m) = ^^(U2)-Ul) + ^U0)) = (C.2) 
^sM = ^ ^5^(3) - 5^(2) + ^s^(l)) = i^t{l + e) , (C.3) 

with t — tan T{m/2N . Since the Fourier transform scales as A"^~^, it depends 
only on the 0(l/iV"^~^) term in Therefore, the higher order transforms 

are, for m < : 

Q^{rn) = -^y|^^(1 + ^')(2 + 3t') , (C.5) 



Using this formalism, it is possible to calculate < Va, d^"'Va > for even n as 

< Va, d~'^Va > = -\ [ [ Va{x)gn{x - x')Va{x')dxdx' 

L Jo Jo 

^ 2N-12N-1 L rL 

= -- V cosjOmCoslOm j / dxMj{x)d:cMi{x')Qn{x - x')dxdx' . 
L ^.^0 Jo Jo 

(C.7) 

This integral contains two types of contributions. Type I arises from the in- 
teraction of distant kinks and type II is a local correction arising from the 
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autocoupling of a given kink and taking into account the discontinuity of 
Qn{x) or its derivatives in a; = 0. At first, wc consider only the type I contri- 
bution for which Qn{^ — x') can be developed as a Taylor series around Xj — xi 
in order to separate the double integration in two independent integrations 
around the kinks 

Qn{x - X') = (^^^-^) + (^^^-^) ((^ - ^j) - - ^l)) 

Now we replace in (C.7) and calculate the local contributions using 

J d^Mjdx = 2r(-l)^' , 
j {x — Xj)dxMjdx = , 
l{x-Xjfd,Mjdx=^{-iy. 



We also expand the trigonometric factor and find that the only non vanishing 
contribution depends on j — I. After relabeling, we obtain 

< >7= E'cOS^(-1>'^„(^ 

that is 

< Va, d.-'^Va >7= -^^nM - ^g^_^{m) + 0(A"-6) (C.8) 

for even n. A similar relation is obtained for odd n : 

A-r2 1 _ 2p2 1 _ 

< d,-^v, - ^-e;„_2(m) + ©(A^"^) . (C.9) 



The local type-II contributions must be examined case by case. For n — 2, 
^2(2^) is continuous in x = but its derivative is not. Near x = 0, we have 



§1 




X 




V6 


L 



The two terms | + are taken into account in type-I contribution. The 
complementary contribution is 

< Va, d-^Va >ii^ J dxM{x)d^M{x')\x - x'\dxdx' . 
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After a bit of algebra, we obtain 

< Va, Va >II^ -JJ 

and finally 

< Va, d-\a >= -y (1 + + ^ ■ (C-10) 

For n — 4:, Qi{x) has a discontinuity on its third derivative in x = 0. It brings 
an 0(1/ A) correction in < Vadx''^Va>, which is of higher order than the terms 
in (C.9). 

For odd orders, including n = 1, the sine factor cancels the type-II contri- 
bution. Notice that for n — 1, Qi{x) is discontinuous in a; = 0. We have 
arbitrarily assumed that ^i(O) = in (C.7) but this is not important since 
only the imaginary part of Qi{m) is used in the sine transform. 

Some other quantities may need to be calculated, of the type < fva, dx~^Va > 
where n is odd and / is a period-A non localized function which is odd over 
the kinks. In this case, the non vanishing contributions are those arising from 
the odd derivatives of Qn- At leading order, we obtain 

or - r 

< fVa, >= -—gn-i(m) j xf(x)d,M{x)dx + ■■■. (C.ll) 

The same expression holds when Va is replaced by Vb. The case n = 1 is special, 
we have 

< fVa, d-^Va >= ^ / f{x)d,M\x)dx . 



D The effect of mean advection 

The case 7 7^ has been studied numerically and near the linear limit Tq (see 
Section 4.1) in Ref. [5]. It is possible to study the stability properties of steady 
solutions for small r in the same way as for 7 = but this is to the price of 
a considerable increase of complexity in the algebra. It is beyond the scope of 
this manuscript to describe the details of these cumbersome calculations. We 
provide here the results without demonstration. 
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When 7 7^ 0, the equihbrium positions of the kinks and antikinks are given by 



, A 

X2p^pA-—, 

A A 
X2p+i = (2p+ 1)- + J, 



where A is related to 7 by 

lere""^ sinh sA 
We define also d — A/A. 



AT 



A 



+ 7. 



(D.l) 



Similarly to Section 4.1, it is convenient to expand the displacements of the 
kinks with respect to the equilibrium in terms of Fourier components. One 
has now to separate the kinks and antikinks as 



N-l 



5x' 



5xi 



2p ^ -rm 
m=Q 
N-l 

H-1 = E ^me' 
m=0 



. (2pfl)m 
N 



By combining these components into 



/ 




we obtain 



^sin^ Or, 



1-C cosh s A - Qr, 




1-CcoshsA + Q^^ 



with 



-sA 



A 
f 

B 



128Ae 
(l-(i2)A^' 

4 

1 - 



4 sin^ ^„ 



d cosh s A — sinh s A 



D 



sAD' 

D = cosh sA — (I sinh sA, 



= (^B + 2C sinh sA + ^C^(cosh 2s A + cos 2^^)^ 



1/2 



(D.2) 
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This result generalizes (18). 



Now, the stabilization effect by friction is still given at first order of the per- 
turbative expansion for small r. We first need to define 

Pm = 2{B cos^ + sin^ 



cos (f)i^^ -^{B cos 0m) , sin = ^ ( sin Om) , 

cos 02m — —(B + C sinh s A sin^ , sin 02m — — —(C sin 2^^ cosh s A) . 

Pm Pm 

The stabilization effect is 

rsin^ 9^ ,2 4 \ 

X ( 1 + d sin 61^ Sin(02m - 01m) - cos Om COs(02m - 01m) " ^ ^^"^ 



(D.3) 

These results have been checked numerically with the same accuracy as those 
presented for 7 = 0. 

In the case of large A and when d is not small, that is when the asymmetry is 
strong, we can neglect all terms 0{e^'^^) in front of terms which are 0(1) or 
larger. Equations (D.2) and (D.3) then simplify considerably. We obtain 

^ 128.3Ae-(A-^) A(l-ci).A)-^0. 

(l + ci)£mA 1 ^ 



a, = -r+-((l-d^)--)(l + d-Ay\ (D.5) 



and 



sKj \ sAj 

As a consequence, the critical value for friction is, at leading order, 

256s3Ae"^(^-^) 



:i + ci)2A 



(D.6) 



Notice that this relation is not valid for small d and does not match (39) for 
d = A = 0. 
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Table 1: Critical values of friction and /3 as a function 
of the number of unstable modes n (or equivalently the 
size of the domain L = 77,(3/2)^/2) and the wavenumber 

of the stationary solution. The columns r^'^ and 
show the analytical predictions given by (39) and (46) 
for 771 = 1. The columns r"^^^ and P^^^ show the critical 
values obtained by numerically solving the perturbation 
problem as indicated in Section 5.1, again for m = 1. The 
columns rj,"™ and show the critical values obtained 
from the direct numerical stability study of the complete 
Cahn-Hilliard equation, as indicated in Section 5.2.2. 



n 


N 


g-sA/2 


„an 
' c 


„pert 
' c 


„num 
' c 


/3r 






7 


2 


410-^ 


6.8478 10-^ 


6.8484 10-^ 


6.84810-^ 


1.0610-4 


2.052210-4 


2.052310-4 




3 


610-3 


2.7469 10-4 


2.7705 10-4 


2.779 10-4 


5.8810-3 


1.632910-2 


1.6385 10-2 




4 


210-2 


5.5655 10-3 


5.914910-3 


6.14710-^ 


5.5810-2 


0.216 


0.222 




5 


510-2 


3.4908 10-2 


3.961410-2 


5.13910-2 


0.260 


1.473 


1.379 


20 


2 


210-10 


6.8832 10-20 


6.8825 10-20 


NA 


3.9310-12 


4.82 10-22 


NA 




3 


410-^ 


4.184610-13 


4.184710-13 


NA 


2.50 10-s 


3.4210-^ 


NA 




4 


210-5 


1.0438 lO-'^ 


1.0439 10-9 


1.00910-9 


2.41 10-6 


3.780 10-6 


3.71610-6 




5 


10-1 


1.1756 IQ-' 


1.175610-^ 


1.1751 IQ-^ 


1.23 10-5 


7.782 10-5 


7.782 10-5 




6 


610-^ 


2.813510-6 


2.814310-6 


2.814410-6 


3.1310-4 


6.85810-4 


6.85810-4 




7 


210-3 


2.7738 10-5 


2.7778 10-5 


2.7795 10-5 


1.4010-3 


3.70610-2 


3.71010-2 




8 


410-3 


1.5681 10-4 


1.575010-4 


1.58010-4 


4.52 10-3 


1.46810-2 


1.47910-2 




9 


710-3 


6.110310-4 


6.1582 10-4 


6.242 10-4 


1.1810-2 


4.676 10-2 


4.82710-2 




10 


10-2 


1.832910-3 


1.841610-3 


1.921 10-3 


2.6810-2 


0.124 


0.138 




11 


210-2 


4.541810-3 


4.4508 10-3 


4.965 10-3 


5.2410-2 


0.287 


0.373 




12 


210-2 


9.7452 10-3 


8.9096 10-3 


1.13510-2 


9.5810-2 


0.551 


1.01 




13 


310-2 


1.8708 10-2 


1.492710-2 


2.392 10-2 


0.164 


0.951 


2.11 
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Table 2: Values of the phase speed of the basic solution 
and the frequency of the m = 1 mode as a function of 
n and N. The columns cl^"^ / (3 and fJ^f^/l3 show the an- 
alytical predictions given by (8) and (37). The columns 
^pert^p and fiX^"^ / j3 show the values obtained by numer- 
ically solving the eigenvalue problem, fi'^ and nX^"^ are 
calculated for m = 1. The columns c^^™/P and 
show the values obtained from the direct numerical sta- 
bility study of the complete Cahn-Hilliard equation near 
the critical value of (3. jj,^'^^ is the frequency of the critical 
mode. 



n 


N 


cT/P 


crV/9 










7 


2 


-17.6492 


-17.6492 


-17.649 


-6.01296 


-6.01298 


-6.0127 




3 


-8.05199 


-8.0485 


-8.0486 


-2.21915 


-2.22148 


-2.222 




4 


-4.61026 


-4.57267 


-4.5730 


-1.16117 


-1.18408 


-1.1876 




5 


-3.11049 


-2.93650 


-2.9368 


-0.63701 


-0.73943 


-0.74685 


20 


2 


-133.043 


-133.043 


-133.04 


-16.4572 


-16.4571 


NA 




3 


-60.8209 


-60.8209 


-60.821 


-5.86179 


-5.86179 


NA 




4 


-35.0024 


-35.0024 


-35.002 


-3.08063 


-3.08063 


-3.0806 




5 


-22.8122 


-22.8122 


-22.812 


-1.92726 


-1.92726 


-1.9273 




6 


-16.0670 


-16.0669 


-16.067 


-1.33119 


-1.33123 


-1.3312 




7 


-11.9316 


-11.9310 


-11.931 


-0.979702 


-0.979894 


-0.97994 




8 


-9.20895 


-9.20716 


-9.2074 


-0.753012 


-0.753716 


-0.75386 




9 


-7.32184 


-7.31636 


-7.3170 


-0.596714 


-0.59863 


-0.59900 




10 


-5.96367 


-5.95021 


-5.9511 


-0.182911 


-0.187182 


-0.18791 




11 


-4.95971 


-4.93139 


-4.9335 


-0.395862 


-0.404176 


-0.40547 




12 


-4.20511 


-4.15177 


-5.1550 


-0.325842 


-0.340601 


-0.34263 




13 


-3.63500 


-3.54226 


-3.5452 


-0.266160 


-0.290808 


-0.2927 
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Table 3: Distribution of solutions as a function of N for 
n = 20, P = and increasing values of r. The statistics is 
calculated out of an ensemble of 100 cases for each value 
of r and A. The cases are obtained by varying the seed 
of the random generator. The fourth and fifth columns 
show the percentages of solutions ending on the stable 
A'^-pair solutions with N in the third column. The first 
column gives the critical linear stability as function of 
A^, provided by the numerical solution to the perturbative 
problem (see Section 3). 



rc 


r 


N 


A ^0.1 


A^l. 


2.81 10"^ 


6 




10"^ 


4 


4% 


19% 






5 


81% 


71% 






6 


15% 


10% 




1.7810-5 


4 


1% 


8% 






5 


57% 


69% 






6 


m 


■2?/X 


2.7810"^ 


7 




3.1610-5 


4 




1% 






5 


44% 


52% 






6 


56% 


41% 




5.62 10-5 


5 


23% 


30% 






6 


74% 


67% 






7 


3% 


3% 




10-4 


5 


9% 


16% 






6 


76% 


75% 






7 


15% 


9% 


1.5710-^ 


8 




1.7810-4 


5 


1% 


2% 






6 


50% 


65% 
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7 


49% 


33% 




3.1610-^ 


6 


IS^X 


31% 






7 


75% 


64% 






8 


7% 


2% 




5.62 10-^ 


6 


1% 


19% 






7 


63% 


62% 






8 


36% 


2% 


6.1610-^ 


9 




10-3 


6 




6% 






7 


29% 


40% 






8 


69% 


54% 






9 


2% 






1.78 10-3 


7 


1% 


20% 






8 


74% 


71% 






9 


25% 


9% 


1.8410-3 


10 




3.1610-3 


7 




3% 






8 


26% 


57% 






9 


68% 


69% 






10 


6% 




4.45 10-3 


11 




5.62 10-3 


8 


2% 


20% 






9 


66% 


70% 






10 


32% 


10% 


8.9010-3 


12 




10-2 


8 




3% 






9 


20% 


51% 






10 


66% 


45% 
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11 


14% 


1% 


1.1910-2 


13 




1.7810-2 


9 


1% 


9% 






10 


33% 


59% 






11 


62% 


32% 






12 


4% 






3.1610-2 


10 


3% 


17% 






11 


48% 


59% 






12 


48% 


23% 






13 


1% 


1% 
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Table 4: Same as Table 4 but for the distribution of solu- 
tions as a function of N for n = 20, r = and increasing 
values of 



beta_c 


beta/3 


N 


A = 0.1 


A=l. 


7.78 10-'^ 


5 




10-4 


3 


12% 


13% 






4 


85% 


87% 






5 


3% 




6.8610-4 


6 




10-3 


4 


17% 


36% 






5 


82% 


63% 






6 


1% 


1% 


3.71 10-2 


7 




10-2 


r> 




2% 






6 


74% 


68% 






7 


26% 


30% 


1.4710-2 


8 


4.68 10-2 


9 




0.1 


7 


2% 


10% 






8 


56% 


78% 






9 


42% 


12% 


0.124 


10 


0.288 


11 


0.551 


12 




1. 


7 




1% 






8 




2% 






9 




17% 






10 


16% 


26% 






11 


65% 


46% 



41 







12 


18% 


8% 






13 
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